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Anomalous Flux Quantization in the Spin-Imbalanced Attractive Hubbard Ring 
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We investigate the one-dimensional Hubbard ring with attractive interaction in the presence 
of imbalanced spin populations by using the exact diagonalization method. The singlet pairing 
correlation function is found to show spatial oscillations with power-law decay as expected 
in the Fulde-Ferrell-Larkin-Ovchinnikov state of a Tomonaga-Luttinger liquid. In the strong 
coupling regime, the system shows an anomalous flux quantization of period h/ie, half of the 
superconducting flux quantum of h/2e, as recently predicted by mean-field analysis, together 
with various flux quanta smaller than h/ie. Notably, the observed flux quanta are determined 
by the difference between the system size Nl and electron number Ne as h/{NL — N^)e. 
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1. Introduction 

Ultracold atomic gases in optical traps liave stimu- 
lated mucli interest in fundamental problems of super- 
fluid states. ^^^-' In particular, two component Fermi 
particle systems with different populations are attract- 
ing much attention due to the possibihty of an exotic 
superfluid phase, which is the so-called Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state.^'^) The FFLO state 
is characterized by the formation of Cooper pairs with fi- 
nite center-of-mass momentum caused by the imbalance 
of the Fermi surfaces of two-component fermions and 
exhibits inhomogeneous superconducting phases with a 
spatially oscillating order parameter. 

Recent theoretical works have revealed that the FFLO 
state accompanies a variety of spatial structures of the 
pairing order parameter, which is dominated by the dif- 
ference between the Fermi wave vectors of the two com- 
ponent fermions. ^^^'*-' On the basis of the bosonization 
approach and conformal field theory, Liischer et al. have 
discussed the correlation exponents of the FFLO state in 
the one-dimensional (ID) attractive Hubbard model and 
have obtained a phase diagram of the electronic states 
on the parameter plane of electron density vs. the spin 
imbalance. ^^^ This diagram indicates that the FFLO re- 
gion expands with increasing strength of the attractive 
interaction. 

More recently, Yoshida and Yanase^^' have pointed 
out that the FFLO state in a mesoscopic ring with 200 
sites exhibits an anomalous flux quantization of period 
/i/4e, which is half of the superconducting flux quantum 
h/2e. They solved the ID attractive Hubbard model on 
the basis of the mean-field approximation in the weak- 
coupling BCS region with the use of the Bogoliubov-de 
Gennes equation.^®' However, the strong quantum fluc- 
tuation effect, which is considered to be crucial for the 
ID mesoscopic ring, together with the strong coupling 
effect in the BCS to Bose-Einstein condensation (BEC) 
crossover region, was not discussed. To achieve a certain 
understanding including such effects beyond the mean- 
field approximation, we believe that nonperturbative and 



reliable approaches are required. 

In the present paper, we investigate the ID Hubbard 
ring with the attractive interaction by using the exact 
diagonalization (ED) method for finite-size systems. To 
confirm the anomalous fiux quantization of period h/Ae, 
we numerically calculate the periodicity of the ground- 
state energy E{^) with respect to the magnetic fiux $ 
without any approximation. Although our calculation is 
restricted to small systems, we expect that the essential 
features of the FFLO state can be well described even 
in finite-size systems, as previously discussed by several 
authors, ^""-'^^) in the case with $ = 0. 

2. Model and Formulation 

We consider the ID Hubbard ring given by the follow- 
ing Hamiltonian: 



H = -tY.i 



gi27r$/JVt^t 



Ci+l,a 



+ /i.c.)-|C/|^n,,tn,4, (1) 



where c] ^ stands for the creation operator for an electron 

with spin a (='f, J.) at site i and Ui^cr — cj^Ci^cr- Here, t 
represents the hopping integral between nearest-neighbor 
sites, and we set t = 1 in this study. $ corresponds to the 
magnetic flux through the ring measured in units of h/e, 
and Nl is the system size. The interaction parameter 
\U\ stands for the strength of the attractive interaction 
on the site. 

We numerically diagonalize the model Hamiltonian 
[eq. (1)] of up to 20 sites using the standard Lanczos 
algorithm. To carry out a systematic calculation, we 
use the periodic boundary condition when A^-|- and N^ 
are odd numbers and the antiperiodic boundary condi- 
tion when they are even, where N-^ and N^ are the total 
numbers of up- and down-spin electrons, respectively.^^^ 
The filling n of electrons is given by n = N^/Nl, where 
Ne{= N^ -\- N^) is the total number of electrons, and 



the spin imbalance is defined by p 



We 



We also 



calculate the correlation function C{r) of singlet super- 



* E-mail address: sano@phen.mie-u.ac.jp 



Kazuhiro Sano and Yoshiaki Ono 




Fig. 1. Singlet pairing correlation functions C'{r) as functions of 
r for several values of p at n=0.4 for \U\ = 10. The solid lines 
represent the results for Ni^ = 20 and the dashed lines represent 
those for N^ = 10. 




Fig. 2. Logarithmic plot of |C(r)| as functions of r for several val- 
ues of p at n=0.4 for \U\ = 10. The solid lines represent the 
results for Nj^ = 20 and the dashed lines represent those for 
Nj^ = 10. The dotted lines are guides for the eyes. 



conducting pairing at the same site as 



3. Pairing Correlation 



(2) 



To confirm the FFLO states and estimate the finite- 
size eflfect in small systems, we examine the singlet pair- 
ing correlation functions C{r) for two different-size sys- 
tems, Nl = 10 and Nl = 20. In Fig. 1, we show C{r) 
as functions of r for p =0 and 0.5 with Nl = 10 and 
for p =0, 0.25, 0.5, and 0.75 with Nl = 20. Compar- 
ing C{r) for p =0 and 0.5 between Nl = 10 and 20, we 
find that the discrepancy between the two sizes is not so 
large. Thus, the finite-size effect is expected to be small 
even in small systems with Nl ^ 10. ^*'' 

We find that C{r) for p = decays monotonically 
as a function of r, while C{r) for p 7^ shows oscil- 
lating behavior. The period of the oscillation seems to 



decrease with increasing p. This result is related to the 
spatial oscillation of the superconducting order parame- 
ter A(r) which has already been pointed out in previous 
works. ^"^®^ The periodicity of A(r) is determined by the 
difference between the Fermi wave numbers of kp-^ and 
kpi and varies as A(r) oc exp(z(5r) for the Fuldc-Fcrrell 
state and as A(r) oc cos{Qr) for the Larkin-Ovchinnikov 



state, where Q 



i-Ft 



«f;| 



mrp 



15) 



In Fig. 1, we 



see that the spatial period of C{r) with n — 0.4 is 20, 
10, and 20/3 for p — 0.25, 0.5, and 0.75, where the cor- 
responding values of Q = nnp are tt/IO, 7r/5, and Stt/IO, 
respectively. Thus, we can conclude that the oscillating 
behavior of C(r) is well accounted for by the FFLO state 
similar to the case of A(r)^"^^' mentioned above. 

Note that, in the purely ID system, the correlation 
function C{r) should show power-law decay for large r, 
as predicted by Tomonaga-Luttinger liquid theory,^^^^^^ 
instead of superconducting long-range order with finite 
A(r) obtained from the mean-field approximation. As 
shown in Fig. 2, the logarithmic plot of |C(r)| indicates 
that the pair correlation function decays as ~ 7.-0.6 ^^^ 
p = 0, and as ^ 7.-11 ^^j. p > 0, in addition to the 
oscillating behavior mentioned before. The present ED 
results with Nl — 10 and Nl ~ 20 seem to be consistent 
with the results obtained from the density-matrix renor- 
malization group"'^^-' for a larger system with Nl = 64.^^^ 
This suggests that the characteristic feature of the FFLO 
state of a Tomonaga-Luttinger liquid is well described 
even in small systems with Nl ^10. 

4. Flux Quantization 

Next, we discuss the flux quantization of the FFLO 
state in the attractive Hubbard ring. Recently, Yoshida 
and Yanase have shown that an anomalous flux quantiza- 
tion with a period of /i/4e, which is half of the supercon- 
ducting flux quantum h/2e, occurs for the FFLO state 
in mesoscopic rings. ^^^ They claimed that the FFLO 
state has a multicomponent order parameter and shows 
rich phases distinguished by the order parameter in the 
presence of the mass flow, which corresponds to the mag- 
netic flux through the ring in our model. When the mass 
flow increases, the FFLO state undergoes a transition 
to another type of FFLO state with a different order 
parameter. Therefore, the ground state of our system 
is expected to alternate between several types of FFLO 
states with increasing of the flux, resulting in anomalous 
flux quantization such as the half period h/Ae. To verify 
this prediction without any approximation, we calculate 
the periodicity of the ground-state energy i<^($) with re- 
spect to the magnetic flux $ by using the ED method. ^^■' 
Although the system is limited to small sizes, this should 
give direct evidence of the anomalous flux quantization 
of the FFLO state, which is well described even in small 
systems with Nl ^ 10 as shown in §3. 

In Figs. 3(a)-3(c), we show the difference in the 
ground-state energy A_E($) = E{<^) — E{0) as a func- 
tion of $ for quarter filling n=0.5 (4 electrons/8 sites) 
at \U\ = 10, 200, and 1000, respectively. In the absence 
of spin imbalance, i.e., p ~ 0, we see that the energy 
levels cross at $ '-^ 0.25 and 0.75 for all values of the at- 
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Fig. 3. Difference in the ground-state energy /\E = E(<^) — E(Q) 
as a function of the magnetic flux $ with p = and 0.5 at 
n=0.5 (4 electrons/8 sites) for \U\ = 10 (a), \U\ = 200 (b), and 
\U\ = 1000 (c). 



tractive interaction \U\, where the usual superconduct- 
ing flux quantization of period h/2e (= 0.5 in the present 
unit of h/e = 1) is observed. On the other hand, in the 
presence of spin-imbalance with p = 0.5, the energy levels 
cross at $ ~ 0.5 for a relatively small value of \U\ = 10, 
as shown in Fig. 3(a), while they cross at $ '^ 0.2, 0.4, 
0.6, and 0.8 for relatively large values of \U\ — 200 and 
1000 as shown in Figs. 3(b) and 3(c), respectively, where 
the period of the ground-state energy with respect to $ 
for p = 0.5 is clearly half of that for p — 0; thus, the 
anomalous flux quantization of period /i/4e takes place. 
The obtained result of the anomalous flux quantization 
seems to agree with that obtained from the mean-field 
approximation.^^) However, the anomalous flux quanti- 
zation is observed exclusively for relatively large \U\ in 
the present study, while it is observed for small |f/| in the 
mean- field study. ^®) This discrepancy will be discussed 
in S5. 
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Fig. 4. Difference in the ground-state energy AE = E{^) — E{0) 
as a function of the magnetic flux $ with p = 0.5 at n=0.5 (4 
electrons/8 sites and 8 electrons/16 sites) for \U\ = 2000. 
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Fig. 5. Difference in the ground-state energy AE = E{^) — E{0) 
as a function of the magnetic flux $ at n = 0.5 (8 electrons/16 
sites) for \U\ = 4000 with p = 0.75 (a), p = 0.5 (b), and p = 0.25 
(c). 



In Fig. 4, we show a comparison of AE{^) between 
two different-size systems, N^ ~ 8 and 16, with p = 0.5 
at n = 0.5 for \U\ = 2000. Here we plot ^ for iV^ = 8 
for comparison with AE for N^ = 16 and find that both 
coincide well with each other for small $. The result sug- 
gests that the size dependence of AE is approximately 
given by 1/Nl for small $. If we assume E{^) to be 
an analytic function with respect to $, E{^) can be ex- 
panded in a power series of $ as 



^($) = E{0) 
1 d^E{^) 



2! 9$2 



$2 



1 d*E{'^) 



4=0 



4! a$4 



$^ + ..., (3) 



b=0 



where the odd powers of $ vanish because of the inver- 
sion symmetry, E{^) = i?(— $). Then, eq. (3) yields for 
small $ 



AE{<i>) ~ 2ttD<^^/Nl, 



(4) 
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Nl=18, p=0.5, |U|=4000 
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Fig. 6. Difference in the ground-state energy A_B = -E($) — £(0) 
as a function of the magnetic fiux $ with p = 0.5 at Nj^ = 18 
and \U\ = 4000 for Ne = 12 (a), A^e = 8 (b), and N4, = 4, (c). 
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$=0 



where D is the Drude weight ;^^'^^^ D = 

and is expected to be independent of Nl for large N^. 
This indicates that the size dependence of Ai<^($) is dom- 
inated by 1/Nl. Therefore, it is considered to be difficult 
to observe the flux quantization for large N^ since the 
value of AE{^) is small for large N^ and finally be- 
comes zero in the limit N^ — > cxd. In fact, Yoshida and 
Yanase^^^ claimed that the anomalous flux quantization 
of h/Ae is observed only in the case with the niesoscopic 
ring and that it disappears in the limit Nl — !> cxd. 

In Fig. 4, we also observe an anomalous flux quanti- 
zation with a period of h/8e, a quarter of the supercon- 
ducting flux quantum h/2e, for the Nl — 16 system, in 
contrast to the Nl = 8 system where a flux quantization 
with a period of h/Ae is observed. It is very curious that 
the two systems exhibit different flux quanta even though 
both have the same n and p, and only the system size Nl 
is different. To clarify this puzzling behavior, we perform 
detailed calculations for various systems. Figure 5 shows 
the $ dependence of /S.E for p = 0.25, 0.5, and 0.75 with 
Nl = 16 and Ne = 8. It clearly indicates that the period 
of these systems is the same regardless of p and is given 
by h/8e. In Fig. 6, we plot the $ dependence of AE for 
several values of N^ with Nl = 18 and p — 0.5, and we 
find that the flux quanta are h/Ae for N^ — 12, h/8e for 
Ne = 8, and h/12e for Ne = 4. These results suggest 
that the flux quanta of the anomalous flux quantization 
exclusively depend on Nl and Ne and are independent of 
p except when p = (usual superconducting state with 
the flux quantum h/2e) and p = 1 (fully spin-polarized 
noninteracting system where the periodicity of A_E($) is 
h/e). 

We perform systematic calculations of the $ depen- 
dence of AE for available systems to determine the flux 
quanta observed for sufficiently large \U\ and obtain the 
phase diagram on the NL-Ne plane shown in Fig. 7, 
where we confirm that the flux quanta are independent 
of p for < p < 1 as mentioned before. Notably, the 
observed flux quanta are determined by the difference 




Fig. 7. Phase diagram of the observed flux quanta on the Nl-Ns 
plane with < p < 1 for sufflcicntly large \U\. 



between the system size Nl and electron number Ne as 
h/{NL — Ne)e as shown in Fig. 7?^^ When we apply 
the result to large systems, we obtain various types of 
anomalous flux quantization including tiny flux quanta 
in the case of large NL—Ne- This is in striking contrast to 
the case with the mean-field approximation, where half 
flux quantum h/Ae is exclusively observed for the meso- 
scopic ring with Nl = 200.^®^ This discrepancy will be 
discussed in §5. 

To qualitatively understand the origin of the anoma- 
lous flux quantization with the notable flux quanta 
h/{NL—Ne)e, let us consider the attractive Hubbard ring 
in the strong coupling limit \U\ ^ 00 for simplicity. In 
this limit, the systems with < p < 1, i.e., A^^^ > A'^j, > 0, 
consist of N^ electron pairs on the same sites (doublons) 
and N-f — N^ unpaired up-spin electrons. When a dou- 
blon moves to the nearest-neighbor site through the hop- 
ping of a down-spin electron of the doublon, an unpaired 
up-spin electron must be at the nearest-neighbor site so 
as not to lose the binding energy —\U\. Therefore, the 
doublons move in the ring accompanied by the unpaired 
up-spin electrons as explicitly shown below. 

Figure 8 shows a schematic diagram of the possible 
hopping process of a doublon to the nearest-neighbor site 
for the simplest system with Ne = 3 {N-f = 2, Ni = 1) 
consisting of a doublon and an unpaired up-spin electron 
in the attractive Hubbard ring with Nl sites in the limit 
\U\ -^ 00. First, we assume the initial state where the 
doublon is at the first site and the up-spin electron is 
at the second site. Next, we move the up-spin electron 
from the second site to the A'^^th site through Nl — 2 
counterclockwise hops. Finally, we move the down-spin 
electron of the doublon from the first site to the A^ith 
site through one clockwise hops to obtain the final state 
where the doublon is at the A'ith site and the up-spin 
electron is at the first site. Then the configuration of the 
final state is rotated in the clockwise direction by one 
site from the initial state. 

The above-mentioned process contains A^^ — 2 counter- 
clockwise electron hops, and one clockwise electron hop, 
corresponding to the Nl — 3 counterclockwise electron 
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Fig. 8. (Color online) Possible hopping process of a doublon to the 
nearest-neighbor site for the simplest case with A^e = 3 {Nf = 2, 



Ni 



1) in the attractive Hubbard ring with N]^ sites in the 



limit \U\ 



hops. Thus, £^(<1>) is considered to exhibit oscillation 
with a period of h/{NL — 3)e as a single electron hops 
under the magnetic flux $ through the ring is respon- 
sible for the oscillation of _E($) with the period of h/e. 
This is in striking contrast to the case of the usual su- 
perconducting state with p = 0, where a doublon (or 
a Cooper pair) hopping process contains two electron 
hops, resulting in the oscillation of E{^) with the period 
of h/2e, that is, the superconducting flux quantum. For 
the general case with N^ electrons where Nf > N^ > 0, a 
similar rotation process contains Nj^ — N-^ counterclock- 
wise electron hops and N_i clockwise electron hops, cor- 
responding to Nl — Ne counterclockwise electron hops. 
Thus, E{^) is considered to show oscillation with a pe- 
riod of h/{NL — Ne)e, which is precisely the anomalous 
flux quantum observed in the present study. 



tive Hubbard model into the repulsive one,^^'^^'^^' both 
systems with the anomalous flux quantization are found 
to be different from each other. ^^^ An interesting future 
problem will be to study the anomalous flux quantiza- 
tion in the attractive Hubbard ring for larger systems 
by using the Bethe ansatz and to elucidate the finite- 
size effect as it might be responsible for the stability of 
the anomalous flux quantization, which is observed ex- 
clusively for the strong coupling regime {\U\ ^ 100) in 
the present ED calculations but for the weak coupling 
regime {\U\ ^ 1.5) in the mean- field approximation.^^' 

Finally, we compare the present results and the previ- 
ous mean-field results concerning the flux quantum. Var- 
ious flux quanta given by h/lN^ — Ne)e are observed 
for the ID attractive Hubbard ring in the present study, 
while only half flux quantum h/Ae was observed in the 
mean-field study^^' where three-dimensionality is implic- 
itly included. Therefore, the present results of the var- 
ious flux quanta except h/Ae might disappear when we 
include the effects of three-dimensionality. In addition, 
although the mean-fleld approximation is appropriate for 
the weak coupling regime, the correlation effects, which 
are neglected in the approximation, become important 
for the strong coupling regime, where the various flux 
quanta are observed in the present study. Thus, the flux 
quanta except /i/4e might be observed with the correla- 
tion effects even in the presence of the three dimension- 
ality in the strong coupling regime. To state our flndings 
more conclusively, it is necessary to investigate the ID 
attractive Hubbard ring including the effects of three- 
dimensionality together with the use of the Bethe ansatz 
available for larger systems. 



5. Summary and Discussion 

We have investigated the ID attractive Hubbard ring 
in the presence of spin imbalance by using the ED 
method, which enables us to obtain reliable results in- 
cluding the strong coupling regime. We have found the 
following. 

(1) The singlet pairing correlation functions show spa- 
tial oscillations with power-law decays as expected in the 
FFLO state of a Tomonaga-Luttinger liquid. 

(2) In the strong coupling regime, the system shows 
an anomalous flux quantization with the flux quantum 
/i/4e, as recently obtained by Yoshida and Yanase^^^' in 
a mean-fleld study, together with various flux quanta 
smaller than h/Ae. 

(3) The observed flux quanta are determined by the 
difference between the system size N^ and electron num- 
ber Ne as h/{NL — Ne)e, independent of the spin imbal- 
ance except when p = and p = 1. 

Although our ED calculations are restricted to small 
systems of up to Nl — 20, the analysis in the strong cou- 
pling limit |C/| — > cxD has confirmed the anomalous flux 
quantization with the notable flux quanta h/^N^ — Ne)e. 
Similar behavior has already been discussed in the repul- 
sive {U > 0) Hubbard ring by using the Bethe ansatz in 
the strong correlation regime with U ^ n\t\, where £'($) 
shows oscillation with a period of /i/7Ve.^^'^^' However, 
by using the well-known transformation of the attrac- 
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